Two-mode optical state truncation and generation of maximally entangled states in 

pumped nonlinear couplers 
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Schemes for optical-state truncation of two cavity modes are analyzed. The systems, referred to as 
the nonlinear quantum scissors devices, comprise two coupled nonlinear oscillators (Kerr nonlinear 
coupler) with one or two of them pumped by external classical fields. It is shown that the quantum 
evolution of the pumped couplers can be closed in a two-qubit Hilbert space spanned by vacuum 
and single-photon states only. Thus, the pumped couplers can behave as a two-qubit system. 
Analysis of time evolution of the quantum entanglement shows that Bell states can be generated. A 
possible implementation of the couplers is suggested in a pumped double-ring cavity with resonantly 
enhanced Kerr nonlinearities in an electromagnetically-induced transparency scheme. The fragility 
of the generated states and their entanglement due to the standard dissipation and phase damping 
are discussed by numerically solving two types of master equations. 



INTRODUCTION 



Methods for preparation and manipulation of nonclas- 
sical states of light have become an important research 
area in quantum optics 0, especially in relation to pos- 
sible optical implementations of quantum computers and 
systems for quantum communication and quantum cryp- 
tography 0. Among various schemes for optical-qubit 
generation, the so-called quantum scissors device of Pegg 
et al. ]jj produces a superposition of vacuum and single- 
photon states, co|0) + ci|l), by optical-state truncation 
of an input single- mode coherent light. The Pegg et al. 
quantum scissors device was studied in numerous papers 



(see, e.g., 
by Babichev et al 




[@>EU)j an d tested experimentally 
111 and Resch et e 
scheme and its generalizations for truncation of an input 
optical state to a superposition of d Fock states (the so- 
called qudits) 0,^3 are based on linear optical elements, 
and thus referred to as the linear quantum scissors de- 
vices. Optical-state 'truncation' can also be achieved in 
systems comprising nonlinear elements (e.g., Kerr me- 
dia) mm, 

and thus will be referred to as the non- 
linear quantum scissors devices. The above-mentioned 
schemes are restricted to the single-mode optical trunca- 
tion. Here, by generalizing our former scheme |l5j . we 
present a realization of nonlinear quantum scissors for 
optical-state truncation of two cavity modes by means of 
a pumped nonlinear coupler. 

Two-mode nonlinear couplers have become, shortly af- 
ter their introduction by Jensen ^(| and Maier one 
of the important topics of photonics due to their wide 
potential a pplic ations and relative simplicity (see, e.g., 
reviews [l8|, |2(|). Among various types of the nonlin- 
ear optical couplers, those based on Kerr effect have at- 
tracted especial interest both in classic al ll^.ll7L fTst ITgj 
and quantum (H El El El El EJUl regimes. 
The Kerr nonlinear couplers can exhibit variations of self- 
trapping, self-modulation and self-switching effects. In 
quantum regime, they can also gen erate sub-Poissonian 
and squeezed light El El til El Possibilities of 
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FIG. 1: A realization of the two-mode nonlinear quantum 
scissors device via a pumped nonlinear coupler in a double- 
ring cavity. Symbols are explained in the text. 



entanglement generation were also studied in nonlinear 
couplers operating by means of the Kerr effect [2ll E3 or 
degenerate parametric down-conversion psj . 



Here, we analyze Kerr nonlinear couplers, which can be 
modelled by systems composed of two quantum nonlinear 
oscillators linearly coupled to each other and placed in- 
side a double-ring cavity. We discuss two schemes based 
on the coupler with an external excitation of a single 
mode and the coupler with two modes pumped. We 
show that the states generated in the excited nonlinear 
couplers under suitable conditions can be limited to a 
superposition of only vacuum and single-photon states, 
c oo|00) +Cqi|01) +Cio|10) +c 11 |11). We compare the pos- 
sibilities of generation of maximally entangled states by 
the couplers excited in single and two modes. We also 
discuss effects of dissipation on the fidelity of truncation 
and suggest a method to achieve strong Kerr interactions 
at low intensities in our system. 
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2. COUPLER PUMPED IN A SINGLE MODE 

We consider a system, referred to as the pumped Kerr 
nonlinear coupler, which consists of two nonlinear oscilla- 
tors, with Kerr nonlinearities Xa and Xb? linearly coupled 
to each other and additionally linearly coupled to an ex- 
ternal classical field. In this section, we assume that the 
field is coupled to one of the oscillators only. Thus, the 
system can be described by the following Hamiltonian 
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and a (b) is the bosonic annihilation operator correspond- 
ing to the mode a (b) of frequency w a Hamiltonian 
fH for H ^l = describes the standard nonlinear coupler 
EElilllllllilllllillil composed of two Kerr non- 
linear oscillators linearly coupled to each other, where 
the parameter e is the strength of this coupling. How- 
ever, our model differs from the standard by inclusion 
of the extra term flgxt, which describes linear coupling 
between the driving single-mode classical field (say, with 

frequency uj^ t ) and the cavity mode a. The parameter 
a describes the strength of this coupling and is propor- 
tional to the classical field amplitude. A possible physical 
realization of the model is presented in figure 1, where an 
external pump, described by flext? is off in the present 



analysis. Kerr media, marked by H n ^ nl and fl no ' nl , are 
linearly coupled as described by grey central region cor- 
responding to flint- 

The evolution of our system in the interaction picture 
can be described by the Schrodinger equation 
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where 



W*)) = ^2 c mn(t)\mn) 



(5) 



(6) 



m,n— 



Thus, the complex probability amplitudes c m „(i) satisfy 
the set of equations of motion: 



«-77 c mn = [Xa™(™ - 1) + Xb™(«- - 1)] C„ 

at 



+ec m _i,„+iV 'm{n + 1) + e*c m +i, n -i\/(m + l)n 
+ac m ^ lin ^/m + a*c m+ i,nVm+ 1. (7) 



A superficial analysis of J7J could lead to a conclusion 
that the evolution of the system pumped by classical ex- 
ternal field cannot be restricted to two lowest photon- 
numer states, but will also include states with a greater 
number of photons. However, by generalizing the method 
of single-mode nonlinear quantum scissors proposed in 
[T^j l (for a review see 1^3, El), we have observed in 
that evolution can be restricted to the only four states: 
1 00) , 1 10} , 1 01) and |11) as a result of degeneracy of Hamil- 

tonians H n ^ nl and H^ nl . By assuming that the couplings 
| a |, and |e| are much smaller than the Kerr nonlineari- 
ties Xa and Xb, we can interpret the evolution between 
the four states as resonant transitions, while the negli- 
gible evolution to other states as out of resonance, anal- 
ogously to the single- mode case 31] . This phenomenon 
can be shown explicitly as follows: under the assump- 
tion of Xa,Xb S> max(|e|, \a\) and short evolution times, 
equation J2J for m,n^0, 1 can be approximated by 

«^c m „ « {x a m(m - 1) + Xb"n{n - 1)] c mn (8) 

which has the simple solution 

c mn {t) w exp{-i [xa«i(m - 1) + Xbn(n - 1)] t}c mn (0). 

(9) 

By setting the initial condition c mn (0) = for m,n ^ 
0, 1, one gets c mn (t) « 0. By contrast, for m,n G {0, 1}, 
the terms proportional to Xa and Xb are vanishing due to 
degeneracy of the Kerr Hamiltonian and so the remain- 
ing terms proportional to e and a are significant. Thus, 
the ideally 'truncated' two-mode state generated in the 
system has the following simple form 

m))cut = c Do (t)|00> + coi(t)|01)+cio(t)|10) + c 1 i(t)|ll), 

where the evolution of c mn , precisely given by J7J), can 
approximately be described by the following equations 
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l dt C ™ 



l dt C11 
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e cio + a en, 



ec i + a c Q. 



(11) 



Hereafter, in equations for the probability amplitudes un- 
der the discussed assumptions, sign '=' should be under- 
stood as Although approximate equations (|TT|l are 
independent of Xa, our derivation clearly shows that the 
Kerr nonlinearity plays a crucial role in the truncation 
process. By assuming that both oscillators are initially 
in vacuum states, \ip(t = 0)) = |00), and parameters a 
and e are real, we find the following solutions of pi(l for 
the time-dependent probability amplitudes: 

c 00 = 77- [(7 - e ) cos ti + (7 + e) cos r 2 ] , 
2 7 
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(f^2 sinri — Q,\ SUIT2) , 
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where Qj = ^2[2a 2 + e 2 + (-l)^-^], 7 = V4a 2 + e 2 , 
and Tj = fljt/2 for j = 1,2. Note that the solution for 
C10 can be written in a more symmetric form since the 
properties hold: (7 + 5)02 = (7— e)Oi = 2^/a 2 (7 2 — e 2 ). 
In a special case of equal couplings a and e, (|12|) simplifies 
to our former solution 15]: 
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where r = at/2. To estimate the quality of the optical- 
state truncation (up to single-photon states) of the gen- 
erated light, we apply fidelity as a measure of discrepancy 
between the ideally truncated two-qubit state p cu t = 
\ip(t)) cv .tcMt{ip(t)\, given by (TJJ, and the actually gener- 
ated output state p(t) = \ip(t)) (ip(t)\ calculated numeri- 
cally from 

\ip(t)) = exp(-itft)|00) (14) 

for a large (practically infinite-dimensional) two-mode 
Hilbcrt space. Specifically, in our numerical analysis, we 
have chosen the dimension equal to 20 for each subspace 
associated with single mode of the field. The fidelity, 
also referred to as the Uhlmann's transition probability 
for mixed states, is defined by (see, e.g., |32^): 

F(p,p cn t) = {Tr[(v/^V^) 1/2 ]} 2 . (15) 

By assuming that one of the states is pure (say |?/>) C ut)> 
then lfH)|l simplifies to F — cu t('/'|p|V ; )cut- Instead fi- 
delity, the Bures distance DB{p\\Pcut) = 2 — 2y/F(p, p C ut) 
is often applied as a measure of discrepancy between the 
states. Note that the Bures distance satisfies the usual 
metric properties including symmetry Db{p\ \f> C ut) = 
-Ds(Pcut| \p), contrary to the quantum Kullback-Leibler 
'distance' (quantum relative entropy) often used in 
quantum- information context j^. General expression 
(|15fl will be applied in our description of the effect of 
damping on the optical-state truncation in section 4. In 
this section, we are focused on pure-state truncation, for 
which equation (|15(l simplifies to the familiar expression 
F = 1 1 (?/;(£) I V^)) cut 1 1 2 , where the ideally truncated state 
|V'(0)cut is given by l|10(l and the actually generated state 
\ip(t)) is calculated numerically from (jl4|l . The fidelity 
for perfect truncation is equal to one. Figure 2 clearly in- 
dicates that the fidelity of truncation using the pumped 
coupler is close to one for short times and the coupling 
strengths much smaller than the nonlinearity parameters 
(|a|,|e| <C Xa,b)- The numerical results shown in figure 
2 confirm the validity of our analytical approach at least 



FIG. 2: Fidelity between the actually generated state IV'(t)) 
and the ideally truncated state |Vw(i)) if the coupling 
strengths are (a) e = a and (b) e = a/10 for the system 
pumped in one mode with /3 = (dashed curves), and that 
pumped in two modes with: /? = a (solid), /? = —a (dot- 
ted), j3 = ia (dot-dashed curves). In figures 2-7, we assume 
that the nonlinearity coefficients are \a — Xb — 10 s rad/sec, 
a — Xa/200, and the coupler is initially in two-mode vacuum 
state 1 00). 



for short evolution times. Thus, we can refer to the sys- 
tem as a kind of nonlinear (as operating by means of Kerr 
nonlinearity) quantum scissors device. 

Solutions for probability amplitudes of the truncated 
states enable a simple calculation of quantum entangle- 
ment, which is one of the most fundamental resources of 
quantum information theory 0. It is well known that 
the entanglement of a bipartite pure state, described by 
a density matrix p — \ip)(ij;\, can be described by the von 
Neumann entropy of either the reduced density matrix 
p a = Tr p or p — Tr a p or, equivalently, by the Shannon 
entropy of the squared Schmidt coefficients pi |32j : 

E(p) = -Tr(p a log 2 p a ) = -Tt(p b log 2 pb) 

N 

= -^2pilag 2 Pi = h(pi,...,p N -i). (16) 

i=i 

For bipartite pure states this measure is often referred 
to as the entropy of entanglement. In a special case of 
two qubits in a pure state, the entropy of entanglement 
E ranges from zero for a separable state to 1 ebit for a 
maximally entangled state, and it is simply given in terms 
of the binary entropy h(p) = — p\og 2 p— (1— p) log 2 (l— p). 
In fact, for a general two-qubit pure state, given by (|1U|) 
with arbitrary amplitudes c mn (m,n = 0, 1), the entropy 
of entanglement given by l|lb|) can simply be calculated 
as 

E(t) = E(\^(t)) cnt ) = E (2|coo(t)cii(t) - c 01 (t)c w (t)\) , 

(17) 
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FIG. 3: Evolution of the entropy of entanglement E of the 
generated states {^(t)) (dots) and the desired truncated states 
IV'cut(i)} (solid curves) by the coupler pumped in a single 
mode with (a) /3 = and in two modes with: (b) f3 — a, 
(c) (3 = —a, (d) f3 — ia, where a = e. 
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FIG. 4: Evolution of the entropy of entanglement E for the 
coupler pumped in (a) single mode (f3 — 0) and (b) two modes 
(/3 = a) using line styles and parameters same as in figure 3 
except e = a/10. On this figure scale, evolutions of E for 
f3 = —a can be described by the same curve as in figure (b), 
while the entropy of entanglement for f3 — ia is negligible as 
less than 5 x 10 -4 . 



where 



(18) 



and h is the binary entropy. If the probability amplitudes 
c mn (t) evolve according to l|12l) then the evolution of the 
entropy of entanglement is given by 

E{t)=E( J2 S{ eco ^3-^)-[e+(-l) J 7 ] 

i=l,2 7 



x cos(ifyi)} sin( 1 ; f73^i)J. (19) 
Solution i|19fl is further simplified by assuming real a — e 
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FIG. 5: Probabilities |fei| 2 for finding the coupler pumped in 
a single mode with e = a in the Bell states: \Bi) and \Bs) 
(solid curves) as well as and \Bi) (dashed curves). 
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FIG. 6: Probabilities for finding the single-mode pumped cou- 
pler with e — a/10 in the Bell states. Same symbols as in 
figure 5. 



then the amplitudes c mn are given by . Thus, we 
obtain 

E(t) = £ f-{[4 + cos(\/5crf)] sin(crf) 

-\/5cos(at)sin(v / 5ai)}). (20) 



Figures 3(a) and 4(a) show the plots of the entropy of en- 
tanglement for the case discussed here, i.e., for the single- 
mode pumped coupler with the coupling parameters |cv| 
and |e| much smaller than the nonlinearities Xa = Xb- 
We see in figure 3(a) that the rapid oscillations in time 
(with a period T\ = 7r/|a|) are modulated by oscillations 
of low frequency (with period Tg = 87r/|a|). As a con- 
sequence, their maxima are of various values but some 
of them approach 1 ebit corresponding to the formation 
of Bell states. To show this explicitly, we represent the 



generated state in the basis 
by the Bell- like states: 



J2j=i b j\ B j) spanned 



= |11)+»100) |fl , = |00)+t|ll 



IS, 



V2 : 
|01) z|10) 
V2 ' 



\Ba 



V2 ■ 
[10> -ij01> 



V2 



(21) 



which differ from the standard Bell states only by the 
phase factor i. Clearly, . \Bj)(Bj\ = 1 and E(\Bj)) = 1 
for i — 1, ...,4, so states (I21|l will shortly be referred to 
as the Bell states. Figures 5 and 6 show the probabili- 
ties for the generation of the Bell states as a function of 
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time [10 5 s] time [10 5 s] 



FIG. 7: Probabilities for finding the coupler pumped in two 
modes with a — f3 and e — a/10 in the Bell states. Same 
symbols as in figure 5. 



time for the single- mode driven couplers, when the ini- 
tial state is the two-mode vacuum state. It is seen that 
states \B{) and \B*i), being superpositions of |00) and 
1 11), are generated with a high accuracy. In detail, the 
maxima of i|20[l occur approximately at times t{m, n) = 
i[(2m-l)2i + (2n- 1)T 2 ] if a = e. Since the frequencies 
of oscillations in (|20|l are incommensurate, by waiting 
long enough (assuming no dissipation), we can achieve 1 
ebit with high precision. For example, in the first four 
periods of T 2 , one could observe generation of entangled 
states approaching Bell state \B\) with E[t(l, 1)] = 0.994 
and E[t(l,2)] = 0.997 ebits, and approaching |£?2) with 
J5[t(2,3)] = 0.992 and E[t{2,A)) = 0.998. Note that only 
the first period is shown in figures 3(a) and 5(a), for 
which the highest entanglement of 0.994 ebits occurs ap- 
proximately at the time i(l, 1) = 4.5 x 10~ 6 . Thus, we 
can effectively treat our system as a source of maximally 
entangled states. On the other hand, the Bell states \B 3 ) 
and I-B4), being superpositions of 1 10) and |01), are not 
generated if the initial state is vacuum and e = a. This 
conclusion can be drawn by observing that the probabil- 
ities for \B%) and I-B4) can reach only 0.8 in figure 5(b). 
However, by relaxing the condition of equal couplings e 
and a, as shown in figure 6(b), Bell states \B 3 ) and |£? 4 ) 
can be generated from vacuum with high precision in the 
dissipation-free system. 

Hitherto, we have assumed that both cavity modes 
were initially in vacuum states. Now, we analyze a more 
general evolution when the cavity modes are initially not 
only in vacuum but also in single-photon Fock states, i.e., 
|^(0)) = (0)> = \kl), where k,l = 0,1. Thus, by as- 
suming as usual that \a\ = \e\ <C Xa,Xb, the evolutions 
of the initial states are found to be 

|V> (01) W)cut = coi|00) + coo|01)+Hii|10) + cio|ll), 
|V> (U V))cut = c 10 |00)+2' 11 |01)+coo|10)+c 01 |ll), 

|^ (11) W)cut = c u |00)+cio|01)+coi|10)+coo|ll) 

(22) 

in terms of the time-dependent amplitudes c mn (r) given 



by lO and 

Coo(t) = cos(v / 5t) cos(t) — sin(\/5"r) sin(r), 
c\i{t) = — i cos(v / 5t) sin(r) — sin(\/5Y) cos(r),(23) 

where, as usual, r = at/2. Please note that Cjj(r) and 
Cjj(r) differ in sign. We find that the generalized ex- 
pression for the entropies of entanglement for the initial 
states \kl) with k, I = 0, 1 reads as 

E^ (t) = S ( -[4 + cos( VEat)] sin(crf) 
V 5 

-(-l) fe "'\/5cos(ai)sm(V5cvt)), (24) 

which implies that 

E V*>)(t) = B(")(t), 

E (01 \t) = E^ 10 \t). (25) 

It is worth noting that our system for the initial states 
1 01} or 1 10) quasi-periodically evolves into the Bell states 
I S3) and I-B4) with high precision but does not evolve 
into \Bi) or \B 2 ) assuming a = e. This is by contrast to 
the evolutions of the initial states |00) (see figure 5) or 
1 11) also for e — a. For brevity, we will not present any 
graphs of the evolutions of the initial states |01), |10), 
and 1 11), which would correspond to figures 2-5 plotted 
for the initial vacua. 

The above solutions for the initial single-photon states 
are included for completeness of our mathematical ap- 
proach to show that, in principle, all Bell states can be 
generated in our system even for e = a. But it should 
be stressed that the system with the initial Fock states is 
much more experimentally challenging than that assum- 
ing initially the vacuum states only. Despite of exper- 
imental difficulty, our system enables generation of the 
one-photon Fock states from vacuum assuming no cou- 
pling e between the modes, which corresponds to having 
two independent pumped cavities with nonlinear Kerr 
media. A possibility of producing single-photon states in 
such systems was demonstrated in [l3| (for a review see 
[3l| and references therein). 

3. COUPLER PUMPED IN TWO MODES 

This section is devoted to the most general scheme 
presented in figure 1, namely that involving two exter- 
nal excitations. We assume here that both modes of the 
coupler are excited by external fields, whereas for the 
case discussed previously we assumed that only one of 
the modes was coupled to the external field. The Hamil- 
tonian describing such a system is of the form 

Hi = + + Hint + H^ t + (26) 

which is the same as that defined by CO~0, except for 
the extra term given by 

H^ t = 0+p*b (27) 
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corresponding to the coupling of the cavity mode b with 
an external driving single-mode classical field (say, with 

frequency Wext)i where the parameter (3 describes the 
strength of this interaction being proportional to the 
classical field amplitude. The evolution of our system, 
described by Hamiltonian (|26() . can be given by the 
Schrodinger equation from which we find the following 
set of equations for the amplitudes c mn (t) of the wave- 
function JfJJl in the interaction picture: 



i^-c m „ = [x a m(m - 1) + Xbn(n - 1)] c n 



+£Cm-i,n+\\/ m{n + 1) + e*c m+ i tn -i\/(m + l)n 

(28) 



+ac m -\, n \fm + a* c m+ i^m + 1 

Vn+l. 



Analogously to the analysis in the former section, we as- 
sume short evolution times as well as the couplings \a\, 
|/3|, and |e| to be much smaller than the Kerr nonlinear- 
ities Xa and Xb- Then, equation (|28|) for m,n/0,l can 
be approximated by © with the solution © , which van- 
ishes for the initial condition c m „(0) = 0. Thus, under 
the above assumptions, the infinite set of equations 128|) 
reduces to the following four equations: 



.da 



(HI 



dt 
.dc m 

l — : — 

dt 
.dc 10 

l — : — 

dt 
.dcu 

% — : — 

dt 



a*c 10 + /3*c i, 
e*ci + a* cn + (3c o, 
ecoi + ac 00 + 13* cn, 
ac i + f3c w . 



(29) 



Assuming that at the time t = 0, both oscillator modes 
are in the vacuum states, i.e., coo = 1 and cqi = cio = 
cn = 0, we can obtain analytical solutions of l|29(l . To 
solve l|29[) we need to find zeros of a fourth-order poly- 
nomial and, hence, the solutions in their general form 
are rather complicated and unreadable. However, if we 
assume that all coupling constants are real and the cou- 
plings with two external fields are of the same strength 
(a = [3) then the solutions become much simpler and 
easier to interpret. Thus, under these assumptions, the 
solutions are found to be: 



1 



coo 



1 



. Xt . 6 , , Xt . 
cos(y) +i-sm{ — ) 

At, 



-iet/2 



C io = ^2^sin(|) e --/ 2 , 
coo - 1, 



(30) 



where we have introduced the effective coupling constant 
A = \/16a 2 + e 2 . As in the case discussed in the previ- 
ous section, the system's dynamics is closed within the 
finite set of n-photon states. Figure 2 shows that for 
the assumed parameters and evolution times shorter than 
5x 10 _6 s, the fidelity between the ideally truncated states 



and the actually generated states by means of the cou- 
pler pumped in two modes deviates from 1 by the values 
less than 0.03 in figure 2(a) or even < 6 x 10~ 4 in figure 
2(b). This again confirms the validity of our analysis and 
justifies referring to this system as a kind of the quan- 
tum scissors device. For the parameters assumed in fig- 
ure 2, truncation with higher fidelity is usually observed 
for the coupler pumped in a single mode rather than in 
two modes. In the latter case, the truncation fidelity de- 
pends on the relative phase between the pumping-field 
couplings a and (3. Note that state iV'cut) for (3 = ±a, 
contrary to (3 = ia, can be directly calculated from so- 
lution (|12|l . For brevity, we have not presented here an 
analogous solution for (3 — ia, but we have used it for 
plotting the corresponding curves in figures 2 and 3(d). 
It is seen, by comparing figures 2(a) and 2(b), that by de- 
creasing e in comparison to ce, the fidelity of truncation 
can be improved for [3 =/= 0. Thus, pumping the coupler in 
two modes can lead sometimes to truncation better than 
that for the system driven in single mode as presented in 
figure 2(b) by dot-dashed curve corresponding to j3 = ia 
and e = a/10. 

The evolution of the pure-state entanglement, given by 
(|16|l . generated in the coupler pumped in two modes can 
be calculated from 1|17JI with the probability amplitudes 
given by (j3t)|l . Thus, we obtain 



E{t)=£[-\l 



A 2 



[16c 



: cos(Ai) + ieAsin(Ai)] 



(31) 

Figures 3(b-d) and 4(b) show the evolution of the en- 
tropy of entanglement measured in ebits as a function of 
time for the system pumped in two modes in comparison 
to the results for the single-mode driven coupler shown 
in figures 3(a) and 4(a). It is seen that the first max- 
imum in figures 3(b-d) is the highest, contrary to the 
case shown in figures 3(a) for the single-mode pumping. 
Nevertheless, the most important fact is that the value 
of the entropy of entanglement E can approach unity to 
a high precision. So, as in the case of the single-mode 
pumping, the two-mode driven system effectively gener- 
ates Bell states. To find which Bell states are generated, 
we can also transform the resulting wavefunction into 
the Bell basis. Thus, figure 7 depicts probabilities for 
the four Bell states. As expected from the form of our 
analytical solutions for c mn (m,n = 0,1), the entangle- 
ment occurs for the states 1 00} and |11), and leads to 
the generation of the states \Bi) and |£?2) (with some 
unimportant global phase factor). Clearly, the highest 
peaks of the entropy E(t) and those of the probabilities 
|fri,2(i)| 2 of the Bell state generation occur at the same 
evolution times, as seen by comparing figures: 3(a) with 
5(a), 4(a) with 6(a), and 4(b) with 7. Analysis of the 
evolutions of the probabilities 1 6x ,2 (*) | 2 and the entropy 
of entanglement E(t) for relatively longer times reveals 
that the oscillations are modulated and some long-time 
oscillations occur in the system, which can be interpreted 
as a result of quantum beats. It is seen from (|31|l that 
two various frequencies appear in our solution and one of 
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them is considerably greater than the other, e.g., the ef- 
fective coupling constant A is VT7 times greater than the 
internal coupling constant of the coupler e = a. We are 
not presenting here dissipation-free evolutions exhibiting 
modulated oscillations at times longer than those in fig- 
ure 4. It would be meaningless since the entanglement 
is lost at such evolution times due to dissipation, which 
inevitably occurs in real physical implementations of the 
coupler as will be discussed in the next section. 

Although our analysis, including all figures, is focused 
on evolution of the initial vacuum states, we present 
shortly some results for other states too. We find that 
the evolution of the initial Fock state |^ feZ )(0)) = \kl) 
with fc, I = 0, 1 is of the form l|22|) but with the proba- 
bility amplitudes c m „(T) given by (|30|) and c mrl (r) equal 
to 



Coo(t) 



1 



-iet/2 



xt 

COS(y) 



^sm(y)]}, 



cn(r) = c 00 (t) - e 1 



(32) 



With the help of these formulae, we can calculate the 
entropies of entanglement explicitly as 

(t) = E 1 1 - A" 2 exp[-z(2|fc - l\ + l)et] 

+ [Wa 2 + e 2 cos{\t)(-l) k - l ie\sm(\t)] ) (33) 



implying the same properties as those given by (|25|) for 
the single-mode excited system. We point out that both 
single- and two-mode pumped couplers for a — e and 
the initial states 1 01} or 1 10) evolve into the Bell states 
I-B3) and I-B4) but do not evolve into \Bi) or |£>2)- This 
is contrary to the evolutions of the initial vacuum states 
(or 1 11)) as, for example, presented in figures 5 and 7. 



4. DISSIPATION 

In a more realistic description both fields a and b lose 
their photons from the cavities. According to the stan- 
dard techniques in theoretical quantum optics, dissipa- 
tion of our system can be modelled by its coupling to 
reservoirs (heat baths) as described by the interaction 
Hamiltonian 

H = Hi + Hi oss , (34) 
Hi oss = f /(a J at)+f 6 /(S,S+)+h.c J (35) 

where 



3 . 



(36) 



are the reservoir operators; c) ' are the boson annihi- 
lation operators of the reservoir oscillators coupled with 
mode a or 6, respectively; g^'^ are the coupling con- 
stants of the interaction with the reservoirs; H\ is given 



either by or l|26(l dependent on the system analyzed. 
We assume two kinds of functions f(a, a)) and f(b, ¥) to 
describe standard damping and dephasing. 

The standard description of a damped system is ob- 
tained for (|35(l with f(a, a') = aJ and f(b,b') — b\ or 
explicitly 



loss 



= f a at + fu + f.st 



(37) 



corresponding to energy transfer between the system and 
reservoirs. It should be stressed that the process de- 
scribed by l|37|) leads to combined effect of amplitude 
and phase damping. The evolution in the Markov ap- 
proximation of the reduced density operator p of the two 
cavity modes after tracing out over the reservoirs can 
be described in the interaction picture by the following 
master equation (see, e.g., [33|) 



dp 
~dt 



-i[Hi_,p] +Ci oss p, 



(38) 



where the LiouviUian 

Aoss/5 = y([ap,a t ] + [a,pa f ]) + -y ([Sp, S+] 

+[S, j 56t])+ 7a n a [[a,p],at]+7 i ,ri 6 [[S,p],6 t ] (39) 

is the usual loss term corresponding to H\ oss , given by 
(|3Tjl : 7fc is the damping rate of the fcth (k = a,b) (ring) 
cavity, and = [exp^Wfc/fcgT) — 1] _1 is the mean num- 
ber of thermal photons at the reservoir temperature T. 
We will analyze both 'noisy'reservoirs (at T > im- 
plying n a ,fii, > 0) and 'quiet' reservoirs (at T 0, so 
n a = fib ~ 0). The latter assumption implies that dif- 
fusion of fluctuations from the reservoirs into the system 
modes is negligible. But still this simplified master equa- 
tion describes the loss of photons from the system modes 
to the reservoirs. 

One can raise some doubts [34| against using Liouvil- 
lian |3IU in a description of lossy anharmonic oscillator 
models given by (J2J. Also the approximation of T = 
in (|39|) is problematic. Nevertheless, master equation 
(|38|l with LiouviUian (|39l) and Hamiltonian Hi set to 

-ffnoni was use d m a number of works both for T > 
(see, e.g.. l35l l36| I37L l3a 1 but also for T = (see, e.g., 
[33I I35I 133 . I39I l40l l41|l. Moreover, the same master 
equation, given by for Hi set to ©, for coupled 
anharmonic oscillators was applied in, e.g.. |37|.l42j. Li- 
ouviUian (|39(l for T — was also used in ( 33] p. 210) to 
describe a model essentially similar to ours comprising a 
nonlinear system, in which two quantized field modes in 
a cavity interact with a classical pump field. The stan- 
dard LiouviUian for T = was also used in to de- 
scribe a system of Kerr nonlinearity, given by H^ nl , and 
a parametric amplifier driven by a pulsed classical field. 
Nevertheless, it should be noted that a realistic master 
equation for Kerr medium [44ll45j is more complicated. 

Phase damping (also referred to as dephasing) can be 
described by l|35|l assuming f(a, dJ) = a^d and f(b, b') — 
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(a) 



time [10"' s] 



(b) 



„0.8 



~0.6 



time [10" 7 s] 

FIG. 8: Effect of the standard damping, described by H38|l and 
139t for quiet reservoirs, on (a) the entanglement of forma- 
tion and (b) fidelity F(p, p C ut), of the states generated in the 
single-mode pumped coupler from the truncated two-qubit 
mixed states for \ a = Xb — 10 8 rad/sec, a = Xa/20, e = a/2, 
and the damping constants f a — ji, equal to (solid), x<j/500 
(dashed), and Xa/200 (dot-dashed curves). Large dots in fig- 
ure (a) correspond to the exact solution. 



b^b, which gives the following loss Hamiltonian 
£io SS = (f a +ft)a+a + (f b + ft)St6. 



(40) 



This interaction can be interpreted as a scattering pro- 
cess, where the number of photons remains unchanged 
contrary to the interaction described by l|37|). Phase 
damping is essential in a fully quantum picture of dis- 
sipation of our system. As a simple generalization of the 
Gardiner- Zoller master equation for a single harmonic os- 
cillator (Hi|, equation (6.1.15)), we describe phase damp- 
ing of our two-mode nonlinear system by master equation 
(|38[l for the Liouvillian 



la 



(2n a + 1) [2a* a pa* a — (a) a) p — p{a) a) 



+y (2n 6 + l)[26 + SpS + S - (b r b) 2 p - p(b*b) 2 }lAl) 

which is significantly different from l|39|) . 

To analyze dissipative evolution, governed by master 
equations l|38|l for Liouvillians (|39|l and <|41|1 . we apply 
standard numerical procedures for solving ordinary dif- 
ferential equations with constant coefficients as an expo- 
nential series. 

The entropy of entanglement E, given by is valid 
for qudits of arbitrary dimension in a pure state, but it 



0.8 



^■0.6 
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UJ0.4 



0.2 



0' 1 — 

0.5 
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1.5 


2 




time [10~ 7 s] 






(b) 

















'0 0.5 1 1.5 2 

time [10" 7 s] 

FIG. 9: Same as in figure 8 but for noisy reservoirs with mean 
number n a — fit — 0.1 of thermal photons. 
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1 

time [10~ 7 s] 
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(b) 
















- 









time [10 s] 

FIG. 10: Effect of the phase damping described by (1381 and 
14111 assuming quiet reservoirs for the same parameters as in 
figure 8. 
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(a) 




time [10" 7 s] 



FIG. 11: Same as in figure 10 but for noisy reservoirs with 
n a = fib = 1, which is 10 times more than in figure 9. 

fails to determine entanglement of a system in a mixed 
state. Thus, for a two-qubit mixed state p, we have to 
apply more general measure, e.g., the Wootters measure 
of entanglement of formation given by |46j| : 

E P (jS)=e(C(p)), (42) 

where £ is given by (|18fl with the argument being the 
concurrence C defined as 

4 

C(p) = max{2maxA 1 - V A l; 0}, (43) 

and Xi are the square roots of the eigenvalues of p{a^ ® 

a^)p*{&^ <g> <7y^), while is the Pauli spin matrix 
of the kth qubit (fc = a, b). It is well-known that the en- 
tanglement of formation goes into the entropy of entan- 
glement for any two-qubit pure states. Examples of evo- 
lution of the entanglement of formation and fidelity for 
the dissipative systems are shown in figures 8-11 both for 
quiet and noisy reservoirs. By analyzing the figures, the 
most important observation is that the generated two- 
qubit entangled states are very fragile to the leakage of 
photons from the cavities in the analyzed couplers. Our 
system is more fragile to losses described by the standard 
master equation, given by H38|) and H39fl . rather than the 
losses due to only phase damping as described by the 
Gardiner-Zoller master equation, given by l|38|) and 1)41(1. 
Inclusion of reservoir noise, at least for the assumed low 
mean numbers n = n a = fit, of thermal photons, does 




FIG. 12: Level structure of atoms in the Schmidt-Imamoglu 
system |4gl 149^1 exhibiting a resonantly enhanced Kerr nonlin- 
earity in mode a. An analogous figure can be drawn for mode 
b. 



not cause a dramatic deterioration of fidelity and entan- 
glement in comparison to the losses caused by coupling 
the system to the zero-temperature reservoirs. Note that 
we have chosen n = 1 in the dephasing model shown in 
figure 11 and much smaller value of n in the standard 
dissipation model in figure 9. The reason is that for the 
latter dissipation model, the number of photons in the 
system can be increased by absorbing thermal photons 
from the reservoirs. If this absorption of thermal photons 
exceeded the loss of the system photons to the reservoirs 
then the generated fields could not be approximated as 
two-qubit states and the Wootters function Ep, given by 
(|42[1 . would fail to be a good entanglement measure. For 
the parameters chosen in figure 9, the generated states 
are well approximated by two-qubit density matrices and 
thus its entanglement of formation is well described by 
(|42[l . By contrast, the number of system photons is not 
affected by the reservoirs in the dephasing model; thus 
(|42H can be used for any number n of thermal photons. 

Finally, it should be noted that fragility of our sys- 
tem to dissipation seems to be a serious drawback from 
an experimental point of view. Nevertheless, a method 
which enables a significant improvement of the entan- 
glement robustness of the generated states, has recently 
been suggested for a similar system |47| . 



5. DISCUSSION AND CONCLUSIONS 

One of the crucial conditions for the successful trun- 
cation and generation of the Bell states in our scheme 
concerns the couplings \a\, \(3\ and |e| to be much smaller 
than the Kerr nonlinearities Xa and Xb- This implies that 
the Kerr interaction should be strong at very low light 
intensities. Thus, it is desirable to discuss an implemen- 
tation, in which such stringent conditions can experimen- 
tally be satisfied. A possible realization can be based 
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on the effect of the electromagnetically-induced trans- 
parency (EIT) or atomic dark resonances as proposed by 
Schmidt and Imamogl u lisl 143 ] (see also |50j ) and ob- 
served experimentally |5lL |52( ■ The Schmidt- Imamoglu 
EIT scheme can be realized in a low density system of 
four-level atoms, which level structure is shown in figure 
12, exhibiting giant resonantly-enhanced Kerr nonlinear- 
ity at very low intensities. The atoms are placed in cavity 
a (and analogously in cavity b) tuned to frequency u> a of 
the mode a resonant with the transition |1) <-*■ |3) and 
detuned by Auj a of the transition |2) <-> |4). The EIT 
effect is created by a classical pumping field of frequency 
Wgxt resonant with the transition |2) <-» |3). By assuming 
|3i3| 2 ^atom/^a < 1 (see |E3)> an ^he atomic levels can 
adiabatically be eliminated, which results in the following 
formula for the Kerr nonlinearity |49| : 



2Xa 



2e V a 



(44) 



present numerical analysis based on the EIT scheme, in 
comparison to 0, E3j more realistic parameters were 
chosen for the coupling constants, Kerr nonlinearities, 
and damping constants. (5) The deterioration of the fi- 
delity of the generated Bell states due to the standard 
dissipation and phase damping was analyzed here within 
two types of master equations both for quiet and noisy 
reservoirs. By contrast, analysis of losses in [l5| was lim- 
ited to the standard master equation and for the quiet 
reservoir only. No effect of dissipation was studied in [43 • 
(6) In the present manuscript, entanglement of formation 
was calculated analytically for dissipation-free systems 
and calculated numerically for dissipative systems. No 
analytical formulae were given in |l5j . while the entan- 
glement of formation was not at all studied in ^tJ- (7) 
The quality of truncation was described here, but it was 
not studied in [TH l47j . Discrepancy between the really 
generated states and the exact truncated states was mea- 
sured here by the fidelity. 



(3) 

where \a is third-order nonlinear susceptibility, gij = 
Hij^Juji/ (2^60^) are the coupling coefficients, tt a is the 
Rabi frequency of the classical driving (and coupling) 
field, Hij is the electric dipole matrix element between 
the states \i) and n a t om is the total number of atoms 
contained in the cavity of volume V a , and eo is the per- 
mittivity of free space. By replacing subscript a by b in 
(|44|l . an analogous expression for \b can be obtained. By 
putting the stringent limit on the required cavity param- 
eters [13, Imamoglu et al. estimated \a ~ 10 8 rad/sec 
|49|. Note that some quantum information applications 
of these giant Kerr nonlinearities have already been stud- 
ied [E EE 13 EE EH- More details about application 
of the Schmidt- Imamoglu scheme of resonantly-enhanced 
Kerr nonlinearity in the analysis of dissipation effects on 
entanglement generation are given by one of us in |57j. 

It is worth stressing the differences between the present 
paper and our former works, (i) In [TH l47j] . only the 
single-mode pumped systems were studied. Here, we an- 
alyze also two-mode pumped systems, which is basically 
a different model, (ii) The model described in (4^ is cru- 
cially different from the one used here as was based on the 
two-mode nonlinear interaction term H- m t — ea 2 H r+h.c. 
In the present work, as well as in [l5|. we apply the two- 
mode linear interaction Hamiltonian H m t — ea^b + h.c, 
which is the same (by neglecting the pump terms -ffext 

and H^ t ) as that used in [13, H3, E3, [13, HE EE E3, S 

The other novelties of the present work in comparison 
to [TEl l47l ] can be summarized as follows: (1) for single- 
mode pumped system, we found a new generalized solu- 
tion l|12fl . which, in a special case of equal couplings a 
and e, simplifies to our former solution obtained in [15j. 
(2) Analytical approximate solutions were found here for 
various initial Fock states |n). In 0,E3> solutions were 
given for initial vacuum states only. (3) Here, we describe 
a possible realization of the model based on the effect 
of the electromagnetically-induced transparency. (4) In 



In conclusion, we have described a realization of the 
generalized two-mode optical-state truncation of two co- 
herent modes via a nonlinear process. Our system is 
a two-mode generalization of the single-mode nonlinear 
quantum scissors device described in 13]. We have de- 
scribed an implementation of the Kerr nonlinear cou- 
plers, where resonantly enhanced nonlinearities can be 
achieved in the Schmidt-Imamoglu EIT scheme. We have 
compared Kerr nonlinear couplers linearly excited in one 
or two modes by external classical fields. We have shown 
under assumption of the coupling strengths to be much 
smaller than the nonlinearity parameters Xa and \b that 
the optical states generated by the couplers are the two- 
qubit truncated states spanned by vacuum and single- 
photon states. Although our approximate solutions of 
the Schrodinger equation are independent of x-constants, 
the Kerr nonlinearity plays a crucial role in the physics as 
we have derived from the complete ^-dependent Hamil- 
tonian. In fact, the Kerr interaction is the mechanism 
in our model, which enables truncation of the gener- 
ated state at some energy level. By contrast, the sys- 
tem without the Kerr nonlinearities and pumped by an 
external field would gain more and more energy. To con- 
firm our predictions, we have compared 'exact' (accurate 
up to double-precision) direct numerical solutions of the 
Schrodinger equation and compared with our approxi- 
mate analytical solutions. The discrepancies between the 
exact and approximate solutions are relatively small as 
shown by fidelities in figures 2-4. We have demonstrated 
that our system initially in vacuum state or single-photon 
Fock states |mn) (m, n = 0, 1) evolves into Bell states. 
We have discussed the fragility of the entanglement of for- 
mation of the generated states due to the standard dissi- 
pation and dephasing in the two distinct master equation 
approaches. 
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